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UNDERGROUND  EXPLOSION  TRIALS  AT  RAUFOSS  1963  : 
BLAST  v'AVE  PROPAGATION  FOLLOWING  A DETONATION 
IN  A TUNNEL  SYSTEM 


SUMMARY 


The  purpose  of  the  present  work  was  to  examine  the  blast 
wave  pressure  to  bo  expected  at  the  entrance  of  a neigh- 
bouring chamber  ir.  an  underground  ammunition  storage  site 
when  al1  the  explosives  in  a chamber  detonate  simultaneously. 

Blast  wave  propagation  and  temperature  distribution  in  a 
simple  tunnel  system  following  a detonation  were  observed 
for  TNT  charges  from  100  to  5400  kg. 

The  blast  wave  pressure  measurements  were  fairly  successful 
whereas  the  temperature  measurements  wexe  not. 

The  peak  pressures  measured  after  firing  a shot  at  the  inner 
end  of  the  tunnel  were  consider  »hly  greater  than  the  corre- 
sponding peak  pressures  measured  after  firing  an  equally 
sized  charge  in  the  chamber.  This  is  due  to  the  branch 
tunnel  having  an  effect  as  a blast  trap. 

A theoretical  model  to  predict  the  expected  peak  pressure 
in  the  transport  tunnel  is  discussed  for  shots  fired  in  the 
chamber.  For  a loading  density  of  200  kg  TNT/m^  an  air 
shock  wave  is  predicted  to  have  an  overpressure  within  a 
range  from  55  bars  to  80  bars.  For  a loading  density  of 
100  kg  TNT/m^  the  overpressure  is  expected  to  be  within 
a range  from  3 5 bars  to  55  bars. 


INTRODUCTION 

Underground  storage  of  ammunition  presents  safety  problems  which  are 
substantially  different  from  those  connected  with  aboveground  storage. 

The  purpose  of  the  present  trials  was  to  examine  several  aspects  of  the 
safety  of  underground  ammunition  storage  with  a particular  view  to 
problems  concerning  the  connected  chamber  type  storage  sites,  which 
are  almost  exclusively  used  in  Norway.  The  trials  were  a continuation 
of  previous  studies  which  have  been  described  in  various  reports 
(Kef  1.  2,  3.  4). 


- . 


The  sp*  cific  objectives  of  the  present  trials  were: 

- to  measure  the  blast  vave  and  the  u mporaturu  di str ; bution  in  a 

tunnel  when  a detonation  occur.,  in  an  aajoining  chamber 

% 

- to  consider  the  blast  wwes  in  view  of  a theory  which  is  described 
elsewhere  in  this  report 

- to  test  a blast  door  and  some  other  equipment  against  the  blast  and 
the  ground  shock  from  a detonation 

- to  measure  peak  pressur  *s  and  pressure- time  histories  in  the  blast 
wav*  outside  the  tunnel 

- to  study  th  amount  and  di  ;tribution  of  debris  ejected  from  the 
tunnel  eXlt 

Sevv  ral  of  the  tests  and  measurements  have  been  reported  elsewhere 
(See  Ref  5,  6,  7,  fi). 

2 BACKGROUND  INFORMATION  /,ND  THEORY 

2 . 1 Expected  general  behaviour 

Some  theoretical  calculations  have  been  made  in  Ref  3 ard  9 concerning 
the  ecuilibrium  pressure  and  temperature  expected  in  a closed  storage 
chamber  when  all  the  explosives  cetonatc  simultaneously.  It  was 
eO*.c,uded  th^t  after  a few  millise*  onds  have  elapsed  to  allow  several 
shock  reflections  in  the  chamber,  an  equilibrium  pressure  of  about 
1070  bars  and  a temperature  of  2100  °K  would  be  attained  for  a loading 
density  of  100  kg  TNT/m\ 

In  Ref  1 it  was  concluded  that  a detonation  in  a chamber  wit.i  a loading 
density  of  100  kg  TNT/m3  will  cause  a shock  wave  of  about  100  bars 
at  the  entrance  to  the  neighbouring  chambers.  If  the  loading  density  is 
increased  to  200  kg  TNT/rn^  it  was  concluded  that  the  a r shock  wave  will 
have  a streneth  of  about  120  bars  t the  entrance  to  the  next  rh^mKor 
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It  was  further  expected  that  the  air  shock  will  be  followed  by  the 
explosive  gases  at  a temperature  between  1000  K a ad  1000  K. 
The  duration  of  the  pressure  profile  would  be  of  th»  order  of  one 
second,  while  the  temperature  duration  would  be  many  seconds 


2 Previous  experiments  in  underground  tunnel  systems  and 

theoretical  considerations 

Several  calculations  and  experiments  concerning  shock  propagation  in 
straight  tunnels  and  shock  tubes  have  been  published  in  the  recent  years. 
However,  most  of  this  work  differs  from  the  present  study  in  several 
ways:  less- than- full- size  smooth  walled  tubes  with  a simple  geometry 
are  considered,  and  the  absence  of  explosive  gases  is  as.  jrru.d.  A theory 
that  would  account  for  the  complex  initial  and  boundary  conditions,  all 
possible  reflections  and  diffractions  of  the  blast  wave,  and  the  energy 
fluctuations  caused  by  tke  shifting  chemical  equilibria,  would  bu  quite 
unmanageable.  Several  simplifying  assumptions  have  thus  to  be  made. 

An  approach  to  the  problem  made  by  D R Curran  and  D Hveding  has  been 
presented  in  Ref  1.  An  exploding  charge  in  a chamber  will  set  up  a very 
complex  system  of  reflected  and  diffracted  shock  waves  which  propagate 
in  different  directions.  The  volume  of  the  chamber  and  the  mass  of  the 
explosive  gases  contained  within  it  is  assumed  to  he  large  compared  with 
the  mass  of  gas  which  passes  through  the  outlet  in  the  first  short  time  porioc* 
necessary  for  some  reflections  to  take  place  in  the  chamber  The  greatest 
pressure  differences  in  the  chamber  are  now  expected  to  be  approximately 
levelled  out.  The  pressure  in  the  chamber  is  then  assumed  to  be  of  the 
same  order  as  the  equilibrium  pressure  which  would  have  been  obtained 
if  the  chamber  was  closed.  This  equilibrium  pressure  and  temperature 
are  calculated  in  Ref  9. 

The  model  is  based  on  the  assumption  that  the  chamber  is  closed  by  an 
imaginary  membran.  After  several  reflections  have  taken  place,  the 
membran  is  allowed  to  burst  caucing  a shock  wave  to  propagate  into  the 
air  of  the  branch  tunnel  (See  Figure  3.1)  and  a rarefaction  wave  to  travel 
back  into  the  explosive  gases.  The  explosive  gas  front  thus  acts  as  a 
piston  supporting  an  air  shock  in  front  of  it.  The  air  shock  overpressure 


** 


•**! 
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i'*  estimated  by  finding  the  inte r Section  of  the  explosive  gas  idiabat  and 
the  Hugoniot  cup*o  of  air. 


In  Ref  11  an  information  summary  of  blast  patterns  in  tunnels  and 
chambers  is  given.  The  data  arc  mainly  based  on  shock  tube  experiments 
in  the  family  of  shock  tubes  at  Ballistic  Research  Laboratories.  In 
Figure  2.1a  curve  from  Ref  11  is  shown.  This  curve  shows  the  incidenr 
3hock  overpressure  versus  the  transmitted  shock  overpressure  for 
a "T"  junction  in  tunnels  with  equal  cros9  section  areas.  The  curve  i a 
supported  by  experiments  up  to  about  10  bars  and  is  accordingly  in 
extrapolation  for  higher  pressures.  This  curve  predicts  an  ippreciablo 
reduction  in  the  blast  wave  pressure. 


In  Ref  12  and  13  shock  propagation  through  channels  and  ducts  is  studied. 
Some  theoretical  calculations  based  on  simplifying  assumptions  and 
experimental  values  are  presented.  The  pressure  are.  however,  small 
and  the  particular  kind  of  junction  of  interest  for  the  present  study  is  net 
included.  Nevertheless,  it  is  clearly  demonstrated  that  a shock  wave 
which  propagates  through  a smooth  channel  with  a cross  section  area 
increase,  decreases  in  strength  and  vice  versa. 


Shock  waves  which  propagate  through  a MTn  junction  will  be  assumed  to 
decrease  in  strength  according  to  Figure  2.1.  If  the  shock  waves  in 
addition  experience  a ross  section  area  increase,  a greater  decrease 
in  the  shock  strength  should  actually  be  the  result. 


Workers  at  Ballistic  Research  Laboratories  (Ref  10)  have  obtained  data 
for  the  propagation  of  shock  profiles  in  a variety  of  tunnel  systems. 
Their  data  were  taken  from  smooth  v/alled  tubes  of  diameters  ranging 
from  a few  inches  up  to  four  feet.  The  attenuation  of  the  peak  press  ir 
in  a long  straight  and  comparatively  smooth  tunned  was  found  to  be  de- 
scribed by  the  formula: 


(p-i) 


(P.-l) exp 


V n 


^ ) KX 
i 


(2.1) 


} 


whore  P is  the  peak  pressure  in  atmospheres  a distance  X from  a point 
where  *he  profile  his  i peak  pressure  ol  P atmospheres  and  a duration 
of  T.  seconds.  D is  the  tunnel  diameter.  K is  given  by 

K-  i . -1— 

* * S U+C  % 

where  S is  the  air  shock  velocity,  U is  the  particle  velocity  behind  the 
air  shock  front  and  C is  the  sound  speed  behind  th*.  air  shock  front. 

The  parameter  Or  is  an  experimentally  determine.*,  parameter  with 
dimensions  of  velocity.  For  smooth  walls  8RL  found  that 
a = (2.93  + 0.262)  m/a,  P is  given  in  atmosp'.  res. 

The  duration  T.  is  defined  as  shown  in  the  sketch  below  and  is  determined 

l 

by  the  slope  of  ’he  pro ssure- time  profile  close  to  the  peak  pressure. 


Tj  -1  Time 


-t/r 

If  the  curve  is  exponential  (P  = P c ' )•  then  it  can  be  shown  that 

T.  is  equal  to  the  time  constant  r. 


In  formula  (2.1)  it  can  bo  shown  that  the  ‘—“term  in  tho  exponential 

^ i KX 

is  dui  to  tho  ov;  rtaking  rarefaction  waves.  The  a yy  term  is  an 
experimental  term  designed  to  account  for  all  irreversible  losc.s, 
which  will  bo  rbitrarily  called  "viscous11  losses. 

It  can  be  seen  that  the  form  of  Equation  (2.1)  is  unfortunate  in  the  case 
of  low  pressure,  since  in  that  case  K — * 0,  and  the  formula  predicts 
no  damping,  whereas  in  fact  we  might  still  expect  viscous  losses.  For 
high  pressures,  it  can  also  be  soon  that  K—  0.  This  is  not  important, 
however,  within  the  pressure  region  which  is  of  interest  here. 

It  should  also  be  mentioned  that  in  Equation  (2.1)  the  term  for  long 
durations  and  for  sufficiently  high  pressures  can  be  small  compared 
with  the  term  y = ^ ^ . This  is  in  accordance  with  the 

expected  behaviour  when  "a  flat"  pressure  pulse  propagates  through 
the  tunnel.  The  attenimtion  is  then  more  likely  due  to  the  "viscous" 
losses  than  to  rarefaction  waves  overtaking  from  behind. 


EXPERIMENTS 


Object 

Tho  specific  objectives  of  the  trials  to  be  described  in  the  present  report 
were  to  measure  the  blast  waves  and  the  tempt:  ratur»  distributions  in 
a tunnel  when  a detonation  occurs  in  an  adjoining  chamber. 

Description  of  experiments 


Experimental  procedure 

The  experimental  procedure  consisted  in  detonating  explosive  charo  s in 
a chamber.  To  achieve  high-pressure  loads  on  a blast  door,  reported  in 
Ref  (6),  charges  wore  also  fired  at  the  end  of  the  tunnel. 


. o 


3.2,2  Tunnel  bystem 


A short  tunnel  - 30  m long  - existed  before  these  trials  were  planned. 

This  tunnel  w is  lengthened  and  two  side  chambers  connected  tc  the 

tunnel  wt-  re  constructed  as  shown  in  Figure  3.1.  The  branch  passage- 

2 

way  had  a cross  section  area  of  2.4  m whereas  the  transport  tunnei 

2 a 

was  5.  S m , T no  volume  of  the  detonation  chamber  was  54  m . 

The  tunnel  system  can  very  roughly  be  considered  as  a half  scale 

model  of  a typical  Norwegian  underground  storage  site. 


The  walls  and  roofs  in  the  two  chambers  were  covered  with  gunite 
(sprayed  concrete), 

The  blast  door  mentioned,  was  installed  at  the  inlet  of  a measurement 
c ham  be  r . 


The  total  di, tonation  chamber  overhead  cover  was  about  13  m.  It  turned 
out,  how.-vc  ;hat  the  soil  layer  in  this  area  was  much  deeper  than 
expected  so  the  thickness  of  rock  w is  only  5-6  rn.  The  strength  of  the 
reck  was  not  very  high  due  to  several  cracks  which  existed  when  the 
trials  started,  and  which  wcr-.  observed  to  be  extended  after  each  shot 
in  the  c ham  be  r . 


3.2.3  Measurement  and  recording  instruments 

The  positions  of  the  measurement  stations  (m  st)  numbered  from  1 to  ^ 
can  be  seen  at  Tigure  3.1  and  Figure  3.2.  Six  measurement  stations 
numbered  from  3 to  3 were  situated  in  the  roof  while  measurement 
station  2 was  situated  in  the  foundation  beside  the  door.  Measurement 
station  1 was  install  *d  in  the*  chamber  behind  the  door. 

'vach  measurement  station  was  equipped  with  a pressure  time  gauge  and 
a thermocouple. 

The  thermocouples  were  Philips  thermocoax  2 ABAC  which  could 
measure  temperatures  up  to  1200  °C  The  response  time  of  a thermo- 
couple is  a function  of  the  particle  velocity  past  the  junction.  The  response 
time  in  milliseconds  for  these  the rmoc ouole s increases  from  ?nn 
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Piezoelectric  pressure  time  transducers  were  used  for  the  pressure 
measurerm  ntr . In  the  tunnel  Kistler  701  H transducers  were  installed 
while  n more  sensitive  Atlantic  Research  Corporation  piezoelectric 
gauge  LC  13  was  mcuntc.l  in  the  chamber  behind  the  door  ( m st  1)  . 

The  Kistler  701  H gauge  had  a sensitivity  of  80  pc/at.  a resonant 
frequency  of  65  kHz.  and  could  measure  pressure  up  to  about  600  bars. 

The  I C 13  gauge  had  a sensitivity  of  500  pc/psi  a.id  could  measure 
pressure  un  to  about  14  bars. 

An  installation  in  the  roof  is  shown  in  Figure  3.3.  The  same  instrument 
holder  wap  used  in  measurement  station  2.  The  cables  from  the  gauges 
ran  through  drill  hole*  to  the  surface.  In  the  measurement  chamber 
behind  the  loor  a different  instrument  hoidcr  bolted  to  a concrete  plat- 
form was  applied.  The  cables  from  measurement  station  Z were  led 
through  a drill  hoY  In  the  coiling  of  the  measurement  chamber  up  fo 
the  surface. 

The  signals  from  the  thermocouples  were  recorded  by  a CEC-galvanometer 
oscillograph.  CFG  7 - 343  galvanometers  were  applied  which  had  a flat 
frequency  range  of  0 - ZOO  Hz  and  a sensitivity  of  0.  165  in/mV. 

The  pressure  gauges  were  connected  to  charge  amplifiers  with  a 
frequency  range  of  0 - 150  kHz. 

The  recording  instruments  were  Ampex  FR  1300  tape  recorders  with 
a frequency  response  from  0 to  20  kHz. 

The  amplifiers  were  placed  above  the  measurement  chamber  on  a 
soft  mat  to  avoid  disturbances  from  the  ground  shock.  The  recording 
instruments  were  placed  a distance  of  about  200  m from  the  amplifiers 
in  a recording  station. 

3,2.4  Experimental  difficulties 

Disturbances  caused  by  the  ground  shock  are  an  inherent  difficulty  in 
underground  trials.  Noice  caused  by  the  ground  shock  was  present  at 
some  pressure  records.  Usually  it  damped  rapidly  out  before  pressures 
were  measured  by  the  gauges. 
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The  cable  leading  from  measurement  station  2 to  its  amplifier  was 
in  most  cases  affected  by  the  ground  shock  during  the  trialo  and  set 
into  mechanical  oscillations . The  pressure  time  history  recorded  at 
this  station  is  therefore  less  reliable.  The  pressure  rise  due  to  iho 
arrival  of  the  blast  wave  should,  however,  be  negligibly  affected. 

Instrument  hollers  as  shown  in  Figure  3.3  were  mounted  in  the  roof 
to  protect  the  gauges  in  the  tunnel.  As  an  undcsired  consequence,  the 
fluid  velocity  past  the  junction  was  too  low  to  provide  sufficient 
response  of  the  thermocouples. 


3 . 3 Experimental  results 

3.3.1  Explosive  charges  fired 

All  charges  consisted  of  cast  TNT  blocks.  An  outline  of  the  shots 
fired  is  given  in  Tabic  3.1  below.  The  positions  of  the  charges, 
the  distance  from  each  charge  center  to  measurement  station  3 
ahead  of  the  door,  and  the  height  above  the  floor  can  also  be  3een 
from  the  table . 


Shot 

Charge 

weight 

kg 

°lacc 

Distance  to 
m st  3 

Distance  from 
charge  center 
to  floor 

i 

100 

Tunnel 

IP..  90  m 

1.40  m 

2 

100 

ii 

ii 

n 

3 

too 

Cham  be  r 

21.15  m 

it 

4 

301 

ii 

21.20  m 

1 . 00  m 

5 

298 

n 

II 

ii 

6 

294 

Tunnel 

1".90  m 

ii 

7 

999 

Chamber 

21.20  m 

0. 50  m 

8 

1003 

Tunne l 

21.85  m 

0.75  m 

9 

1001 

ii 

17.68  m 

0.  50  m 

10 

5426 

Che  mber 

21 .90  m 

*• 

Tabic  3 . 1 Shots  fired 


f i q f v <• 
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3.3.2  Calculated  pressure  based  on  the  arrival  times 

The  time  for  the  blast  wave  to  arrive  at  the  respective  measurement 
stations  was  measured.  The  number  of  measurements  for  each  shot 
varie  d from  4 to  6.  The  data  are  presented  in  Table  3.2.  A smooth 
curve  was  then  fitted  to  the  arrival  time  distance  data  by  the  method 
of  least  squares.  A discussion  of  choice  of  curves  and  method  of 
calculation  is  given  in  Ref  (14).  The  shock  arrival  time  versus  distance 
from  the  charge  obtained  by  curve  fitting  to  the  measured  data  can  bo 
seen  in  Figure  3.4  and  Table  3.3.  It  should  be  mentioned  that  the 
time  of  trigger  was  lost  for  shot  4,  5 and  7.  For  these  shots  the 
arrival  time  at  measurement  "tation  3 was  arbitrarily  chosen  as 
20  ms.  The  curves  wore  diffo rentiated  and  the  blast  wave  velocity 
obtained.  The  v-locity  and  the  appropriate  Rankine-Hugoniot  relation 
determine  the  blast  wave  pressure  along  the  tunnel. 

Pressure-distance  curves  could  either  be  based  on  the  distance  from 
the  charge  or  from  a fixed  point  in  the  tunnel,  for  instance  m st  3 . 

Both  type  of  curves,  where  pressure  means  the  absolute  pressure 
(chock  overpressure  plus  atmospheric  pressure),  are  presented  in 
Figure  3.5  - 3 . 8. 

3.3.3  Measured  pressure 

The  pressure  records  from  measurement  stations  2-8  are  presented 
in  Figures  3.9  - 3.23.  The  pressures  calculated  from  arrival  time 
data  arc  also  shown  on  the  records.  It  should  he  noted  that  the  difference 
in  the  calculated  pressure  between  measurement  station  2 and  3 due  to 
their  position  along  the  tunnel  is  small,  and  the  calculated  pressures  f^r 
measurement  station  3 are  therefore  also  indicated  at  the  records  from 
measurement  station  2.  The  pressure  records  presented,  are  traced 
from  the  original  records  and  some  loss  in  quality  may  therefore  be 
expected. 

The  quality  of  the  measurements  was  not  constant  for  all  shots.  Drift 
in  the  zero  adjustment  can  be  seen  on  records  from  the  first  five  shots. 
This  difficulty  was  avoided  for  the  later,  shots.  The  precision  of  the 
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measurement  chain  was  estimated  io  about  10  %.  The  initial  values  of 
the  pressure  curves  are  thus  expected  to  be  of  a similar  precision. 

The  precision  of  the  pressure-time  history  is  also  expected  to  be  quit< 
good  when  the  zero  adjustment  of  the  amplifiers  is  proper.  The  pre- 
cision at  later  stages  of  the  pressure  records  may,  however,  be 
slightly  reduced  due  to  heat  exchange  from  the  explosive  gases  to  thj 
gauges.  In  view  of  the  uncertainty  in  the  interpretation  of  the  pressure 
curves,  the  precision  is  considered  to  be  sufficient. 

The  pressure-time  history  recorded  from  measurement  station  2 wat 
generally  of  poorer  quality  than  the  pressure  measurements  ai  other 
stations.  This  was  due  to  vibrations  in  the  cable  leading  to  the  gauge. 
Several  attempts  to  prevent  this  were  not  successful.  The  records 
from  measurement  station  2 are  presented,  however,  because  the  rapid 
rise  in  pressure  duo  to  the  arrival  of  the  shock  front  is  assumed  to  be 
relatively  unaffected  by  the  induced  disturbance  in  the  cable. 

The  pressure  records  from  shot  10  differ  from  records  from  the  previous 
shots,  "Hie  overhead  cover  of  the  detonation  chamber  was  not  of  good 
quality  and  sufficient  thick  to  contain  the  effects  of  this  shat,  and  a 
crater  was  formed.  As  a consequence  the  pressure-time  history  in 
the  tunnel  was  highly  affected  by  a decrease  in  the  duration.  The 
"piston"  * ffect  supporting  the  shock  propagation  will  be  weaker  than  in 
a contained  detonation.  Accordingly  the  attenuation  of  the  blast  wave 
will  increase  more  with  distance  along  the  tunnel  in  this  case.  Another 
effect  from  the  blow-through  was  that  the  amnlitiers  and  cables  at  th«. 
surface  were  exposed  to  considerable  disturbances  by  moving  rock  Jrciv, 
the  crater  which  broke  cables  and  hurried  the  amplifiers.  Nevertheless, 
as  a crater  formation  takes  some  time,  the  pressure-time  history 
recorded  at  m ?t  2 and  3 within  the  first  10  milliseconds  is  likely  to 
be  unaffected  by  the  formation  of  a crater. 

3,3.4  Discussion  of  pressure  records 

Pressure  records  from  pressure  measurements  in  the  open  are  often 
approximated  by  exponential  curves.  Pressure  curves  from  measure- 
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ments  o£  the  shock  wave  propagation  in  smooth  pipes  and  in  shock  tubes 
may  also  frequently  exhibit  exponential  shape.  The  pressure  records 
presented  here  differ  in  some  respects  from  exponential  shape.  They 
start  with  the  ordinary  steep  rise  at  the  initial  stage  due  to  the  shock 
front,  immediately  followed  by  some  high  and  narrow  spikes  behind 
which  broader  spikes  take  place  before  the  curve  starts  to  smooth  out 
as  the  overpressure  decreases. 

In  previous  experiments  in  underground  tunnel  systems  (P^cf  l)  the 
pressure-time  elapse  was  observed  to  be  somewhat  different  from 
what  is  observed  here.  A ^ow  pressure  precursor  followed  by  several 
peaks  of  higher  pressure  was  measured,  and  the  shock  velocity  usually 
was  in  agreement  with  the  low  pressure  precursor.  It  was  thought  that 
the  low  pressure  precursor  was  due  to  the  spherically  expanding  part 
of  the  original  blast  wave,  which  at  the  moment  was  not  overtaken  by 
the  higher  pressure  waves  which  were  built  up  from  behind  through  the 
process  of  multiple  reflections  from  the  walls. 

According  to  this  conception  it  follows  from  the  initial  parts  of  the 
pressure  curves,  Figures  3.9  - 3.23,  that  the  measurement  stations 
have  been  at  sufficient  distance  from  the  charges  to  ]ct  the  higher 
pressure  waves  due  to  multiple  reflections,  overtake  the  expanding 
part  of  the  original  blast  wave. 

It  is  then  likely  that  the  blast  wave  following  a detonation  in  the  tunnel 
system  consists  of  a relatively  plane  shock  wave  moving  along  the  tunnel 
axis,  superposed  by  a great  number  of  waves  generated  by  the  process 
of  multiple  reflections  and  diffractions  with  the  irregular  walls  in  the 
tunnel.  The  calculated  pressure  based  on  the  arrival  times  should  then 
approximate  the  plane  shock  wave  component, which  may  frequently 
have  a lower  value  than  the  total  peak  overpressure  felt  by  the  gauge. 

Local  geometrical  irregularities  may  also  be  responsible  for  some  of 
the  anomalies  of  the  pressure  records.  As  an  example  consider 
measurement  station  3.  Figures  3.2  and  3.24  show  the  position  of  the 
station.  During  the  excavation  of  the  measurement  chamber  a mass  of 
heavy  rock  in  the  roof  ahead  of  the  chamber  inlet  fell  down.  To  mount 
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the  station  it  the  level  of  the  roof  profile  a protruding  concrete  foundation 
was  made.  No  calculation  can  be  made  to  show  the  influence  of  this  to 
the  records.  3ut  the  records  from  this  station  exhibit  stronger  oscil- 
lations at  the  initial  stages  thar  records  from  other  stations  in  the  roof. 
The  records  from  this  station  also  exhibit  higher  peak  pressures  than 
records  from  measurement  station  2. 

The  measured  peak  pressure  is  according  to  the  previous  discussion 
expected  to  be  equal  to  or  greater  than  the  calculated  pressure  (derived 
from  arrival  time  data).  The  maximum  side-on  pressure  in  the  blast 
wave  should  generally  be  bounded  by  the  calculated  pressure  and  the  peak 
pressure  felt  by  the  gauge.  tvhen  an  appreciable  difference  exists  between 
the  measured  peak  pressure  and  the  calculated  pressure,  the  sidc-on 
pressure  in  the  blast  wave  may  be  approximated  by  smoothening  the 
pressure  curves  to  for  example  an  exponential  shape  by  conserving  the 
“re a under  the  .urve. 

Accordingly  the  pressure  curves  presented  here  will  bo  interpreted 
in  the  following  ways: 

1)  The  maximum  pressure  c in  be  taken  directly  from  each  curve 
when  some  care  is  taken  to  avoid  spikes  which  do  not  contain 

an  appreciable  area.  Such  spikes  could  possibly  be  due  to  local 
disturbances  near  the  gauge  and  to  noice  in  the  measurement  chain. 

2)  The  curves  can  be  :moothened  to  an  exponential  shape  so  that  the 
areas  under  the  curves  are  approximately  conserved. 

The  blast  waves  from  detonations  in  the  chamber  were  of  main  importance 
for  the  purpose  of  this  report.  Some  shots  were,  however,  fired  at  the 
end  of  the  tunnel  particularly  to  expose  a blast  door  to  high  pressure  blast 
wave  loads  (Ref  6).  These  shots  are  also  included  in  Tabic  3.1.  A look 
at  the  figures  containing  the  calcu'atcd  pressure  (Figures  3.5  - 3.8)  and 
the  pressure  records  (Figures  3.9  - 3.23)  demonstrates  that  the  measured 
pressures  from  shots  of  equal  charge  size  fired  in  the  chamber  and  in 
the  tunnel  are  different. 
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At  m st  3 the  peak  pressure  of  the  blast  wave  following  a charge  fired 
at  the  end  of  the  tunnel  is  considerably  greater  than  the  blast  wave 
following  a similar  charge  fired  in  the  chamber.  This  difference  se^ms 
to  increase  with  the  charge  9ize. 

Blast  waves  from  charges  fired  at  the  end  of  the  tunnel  tend  to  be  of 
shorter  duration  compared  with  similar  shots  in  the  chamber  as  the 
pressure  records  exhibit  a more  rapid  decay  with  time.  Pressure- 
time  histories  from  tunnel  detonations  also  tend  to  have  smoother 
pressure  profiles  than  the  pressure  profiles  following  chamber  deto- 
nations. This  is  frequently  the  case  for  pressure  records  from  m st  3. 

Interesting  differences  are  observed  when  we  compare  the  pressure  data 
from  shots  8 and  9.  The  charge  fired  as  shot  8 was  arranged  in  a stack 
close  to  the  end  wall  of  the  tunnel.  Figure  3.25.  This  shot  resulted  in 
several  cracks  in  the  tunnel  walls  and  the  roof  near  to  the  charge.  For 
safety  reasons  the  following  shot  was  fired  about  4.2  m from  the  end  of 
the  tunnel  (Tabic  3.1). 

The  difference  in  peak  pressure  is  considerable  in  this  case  as  dis- 
played for  example  on  Figures  3.6  and  3.8.  The  attenuatior  differs 
also  between  the  two  shots  as  the  blast  wave  propagates  through  the 
tunnel.  The  peak  pressure  in  the  blast  wave  following  shot  8 is  much 
greater  than  the  corresponding  pressure  following  shot  9.  The  peak 
prcSHure  in  the  blast  wave  following  shot  8 attenuates  stronger  than 
what  is  the  case  after  firing  shot  9,  ar.d  it  decreases  below  the  corre- 
sponding peak  pressure  in  the  blast  wave  of  the  following  shot  near 
the  end  of  the  tunnel. 

3.3.5  Pressure  duration 

Both  the  definition  and  measurement  of  the  pressure  duration  introduce 
difficulties.  Some  of  the  pressure*  time  profiles  measured  here  can 
approximately  be  considered  to  have  triangular  shape.  Frequently, 
however,  it  is  better  to  approximate  the  pressure  curves  by  an  expo- 
nential curve  at  the  initial  stage.  It  follows  that  the  term  duration 
would  mean  something  quite  different  in  these  two  cases. 
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Consequently,  the  term  duration  must  be  defined  to  be  meaningful.  But 
the  duration  according  to  a definition  can  not  be  found  independently  of 
how  the  pressure  curve  is  road  md  the  side-on  pressure  assigned  to  the 
blast  wave  is  derived. 

Accordingly  any  sot  of  values  of  a duration  will  contain  great  uncertainties. 
A set  of  such  data  will  further  not  provide  any  measure  of  time  in  which 
appreciable  pressure  can  be  measured.  The  term  duration  is  not 
particularly  useful,  «\r.d  it  is  thus  hard  to  find  support  to  present  a set 
of  such  data  from  the  pressure  records. 

3.3.6  Temperature  measurements 

As  mentioned  elsewhere,  the  thermogauges  at  each  measurement  station 
in  the  roof  were  locatv  d in  an  installation  which  is  shov'n  in  Figure  3.3. 

The  installation  provided  necessary  protection  of  the  gauge  from  the 
blast  waves.  As  an  uncle  sired  consequence,  the  fluid  velocity  past 
the  junction  was  too  low  to  provide  sufficient  sensitivity  of  the  gauge. 
Temperature  deflection  was  recorded  it  the  galvanometer  oscillograph, 
but  the  maximum  temperature  measured  was  far  too  low,  and  no  cor- 
rection can  be  given.  The  temperature  mease remonts  were  thus  un- 
successful. The  temperature  tecords  following  shot  7 are  as  an 
example  snown  in  Figure  3.27. 


4 COMPA RISQN^O F THEORY  AND  EXPERIMENTS 

4 . 1 Tentcativc  correction  of  thi  BRL  data 

The  exponent  in  Ho  2.1  contains  the  factor  a = (2.93  P^  + 0.262)  m/s 
which  wa3  experimentally  determined  in  smooth  walled  tubes. 

The  experimental  constants  in  the  expression  for  a , Eq  2.1  will  here 
be  tentatively  adjusted  according  t<»  the  present  condition  with  rough 
walls  etc.  The  value  of  c/ obtained  that  way  will  then  bo  applied  to 
the  model  previously  described  to  predict  the  pressure  in  similar 
tunnel  systemr. , 
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The  experiments  of  the  present  trials  differed  from  those  of  BRL  in 
two  important  respects: 

1)  The  presence  of  hot  explosive  Rases  in  the  tunnel 

2)  The  roughness  of  the  tunnel  walls 

As  stated  earlier  the  term  containing  or  was  designed  to  take  account 
of  all  the  irreversible  losses.  Accordingly  it  is  expected  that  the 
experimentally  determined  constants  contained  in  a would  have  been 
different  if  the  walls  in  the  shock  tubes  had  been  rough.  The  constant 
multiplied  by  P9  in  the  expression  for  a contributes  more  to  a than 
the  additive  constant  does  even  for  moderate  pressures,  so  the  latter 
will  be  neglected  here. 

It  was  previously  discussed  why  the  pressure  duration  was  difficult  to 
account  for.  A definition  of  the  pressure  duration,  as  the  time  for  the 
pressure  to  drop  to  l/e  of  the  maximum  value,  applied  to  these  curves 
might  result  in  values  in  the  range  of  20  - 100  ms.  An  inspection  of 
the  exponent  ir.  Eq  2, 1 shows  that  the  term  containing  a may  particu- 
larly for  blast  waves  of  higher  pressure  and  long  duration  in  a tunnel 
with  rough  walls,  dominate  the  term  containing  l/T. . It  is  there- 
fore worthwhile  to  see  if  the  pressure  data  can  provide  support  to 
this  assumption. 

Pj  and  P^  represent  the  calculated  pressure  at  measurements  station 
3 and  7 respectivclv,  in  Table  4.1  below.  The  j.^ak  pressures  from 
the  records  according  to  interpretation  2 might  in  principle  as  well 
have  been  chosen.  The  peak  pressure  is,  however,  not  as  well  defined 
as  the  calculated  pressure.  This  is  unfortunate  as  the  term 
( ^ + l/Tj)  appears  to  be  sensitive,  when  it  is  calculated,  to  uncer- 
tainties in  the  pressure  data.  and  in  the  table  are  the  distances 
respectively,  from  the  charge  to  measurement  station  3 and  7. 
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Q 1 

The  calculation  of  ( - + - ) i3  based  on  Eq  2. 1 and  data  presented  in 

i 1 c/ 

the  table.  when  it  is  assumed  that  j,.  can  be  ignored  relative  to 

the  ratio  •=>  can  be  calculated  and  incorporated  in  the  table. 


An  inspection  of  Table  4.1  shows  that  the  ratio  ^ is  in  the  range  cf 

20  - 30  for  the  300  kg,  1000  kg  and  5400  kg  shots  which  provides  some 

support  to  the  assumption.  If  durations  in  the  range  of  3 - 35  rr.3  arc 

assumed  with  shorter  duration  to  higher  pressures,  a significant 

decrease  in  the  values  of  ^ , presented  in  the  table  is  attained. 

*1 

A greater  number  of  measurements  and  shots  than  attained  in  these 
trials  is  necessary,  however,  to  verify  these  ideas.  Nevertheless, 

1 fy 

if  we  ignore  the  term  — and  apply  as  an  effective  value  — = 20 

1 i 

for  shots  fired  in  the  chamber,  a conservative  estimate  of  the  atte- 
nuation would  be  the  result  for  higher  pressures. 


4 . 2 Data  obtained  from  the  theoretical  model 


The  detonation  chamber  is  assumed  to  be  closed  and  its  initial  pressure 
depending  on  the  charge  density  is  calculated  by  the  methods  of  Ref  (9). 
The  results  are  presented  in  Table  4.2  (page  28)  for  actual  shots  and 
charge  densities. 


The  air  shock  which  propagates  into  the  branch  tunnel  is  then  calcu- 
lated by  fimiing  the  intersection  between  the  gas  adiabat  for  the  explo- 
sive gas  and  the  Hugomot  curve  of  air  for  each  shot. 


This  is  done  for  y = 1.3  and  y =1.5.  The  first  value  of  y is  to  be 
applied  for  moderate  pressures  while  y = 1.5  is  supposed  to  be  an 
effective  value  of  y to  account  for  real  gas  effects  at  higher  pressures 
and  temperatures.  The  results  are  presented  in  the  table. 


The  shocks  which  propagate  down  the  5 m long  branch  tunnel  are 
assumed  to  attenuate  according  to  the  DRL  formula  with  an  effective 
value  of  ^ ~ 20  thus  providing  a conservative  estimate  of  the  atte- 
nuation. Values  for  the  blast  wave  peak  pressure  at  the  outlet  of  the 
branch  tunnel  estimated  by  this  method  arc  also  presented  in  Table  4.2, 


These  data  arc  finally  applied  to  the  extrapolated  data  of  Ref  (11)  which 
are  presented  in  F.gure  2.1.  The  results  can  be  seen  in  the  table  for 
both  values  of  y together  with  the  measured  peak  pressure  from  measure 
mont  station  3 which  was  the  station  situated  nearest  the  junction. 


Discussion 


The  theoretical  model  presupposes  that  the  chamber  is  closed  by  an 
imaginary  membran  until  equilibrium  pressure  is  attained.  For  a 
loading  density  of  100  kg  TNT/m^  the  equilibrium  pressure  is  esti- 
mated to  1070  bars.  As  a comparison  we  will  consider  the  blast 
wave  from  a centrally  located  detonating  charge  in  a spherical  chamber 
as  the  blast  wave  reaches  the  chamber  wall.  For  the  same  loading 
density  the  leading  blast  pressure  at  this  distance  i3  250  bars.  This 
shows  that  the  wave  front  has  a much  lower  pressure  than  the  gases 
behind  which  have  a maximum  pressure  probably  well  above  1070  bars. 
It  is  then  also  clear  that  a classical  shock  front  with  both  a uniform 
decaying  pro  ssure -time  history  and  a uniform  decaying  pressure- 
distance  distribution  is  not  yet  formed  at  this  stage.  The  comparison 
indicates,  however,  that  the  model  is  a fair  approximation  of  the 
actual  phenomena. 


The  use  of  the  extrapolated  curve  in  Figure  2,1  represents  some 
uncertainty.  More  important  is  probably  the  fact  that  the  cross 
section  area  of  the  transport  tunnel  is  twice  the  size  of  the  branch 
tunnel.  The  attenuation  of  the  blast  wave  in  the  MTM  junction  seems 
therefore  to  be  underestimated. 

Some  attenuation  of  the  blast  wave  will  also  occur  during  its  travel 
from  the  function  to  measurement  station  3,  which  represents  a 
distance  of  11  m. 

A disturbing  factor  is  the  reflection  of  the  blast  wave  £t  the  end  of 
the  transport  tunnel,  which  would  increase  the  peak  pressure  of  the 
actual  wave  compared  to  that  of  the  model.  It  does  not  seem  likely, 
however,  that  this  effect  is  appreciable  at  such  a short  range  as  the 
distance  from  the  junction  to  measurement  station  3. 

Table  4.2  shows  a fair  agreement  between  the  measured  peak 
pressures  according  to  interpretation  1 and  those  obtained  by  the 
model  using  y s 1.5.  It  may  therefore  be  expected  that  the  model 
could  also  be  used  with  some  confidence  at  higher  loading  densities 
where  experimental  data  are  lacking.  This  has  been  done  in  Table  4.2 
for  the  loading  densities  200  and  27C  kg/m^. 

It  is  obvious,  however,  that  considerable  uncertainties  are  involved 
in  the  methods  used  for  the  theoretical  calculation.  Also  the  experi- 
mental data  should  be  regarded  with  some  suspicion,  particularly 
at  the  shortest  distances,  where  unpredictable  reflections  and  re- 
fractions may  cause  a very  irregular  propagation  of  the  blast  wave 
with  great  local  variations  of  the  blast  pressure.  Previous  trials 
(Ref  l)  indicated  higher  pressures  than  those  reported  here.  Even 
if  the  present  results  are  thought  to  be  more  reliable  due  to  an  im- 
proved measurement  technique,  it  is  nevertheless  felt  that  some 
consideration  should  be  given  to  the  previous  results.  It  is  there- 
fore suggested  that  the  present  values  of  the  sidc-on  pressures  should 
be  multiplied  by  a factor  of  1.5  to  obtain  an  estimate  of  the  upper 
limit  of  the  shock  strength  for  the  various  loading  densities. 
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The  model  exhibits  some  properties  which  should  bo  noted. 

The  blast  wave  attenuates  while  it  propagates  through  the  branch  tunnel. 
This  can  be  soon  from  Table  4.2  for  both  values  ofy.  Table  4.2  also 
reveals  that  the  attenuation  increases  with  pressure.  Eq.  2.1  shows 
that  the  pressure  attenuation  also  increases  with  increasing  length  and 
reciprocal  diameter  of  the  branch  tunnel.  Rough  walls  contribute  in 
the  same  direction.  Tie  numbers  in  the  paranthesis  indicate  as  an 
example  what  to  be  expected  when  the  cross-section  area  of  the  branch 
tunnel  is  5 m^. 

Further  attenuation  takes  place  as  the  blast  wave  propagates  through 
the  junction  connecting  the  branch  tunnel  and  the  transport  tunnel. 

The  data  presented  in  Figure  2.1  are  given  for  tunnels  of  equal  cross- 
section  areas.  An  increase  in  the  ratio  between  the  cross-section 
area  of  the  transport  and  the  branch  tunnel  would  provide  a greater 
pressure  decrease  than  indicated  by  the  data  of  Figure  2.1.  For  the  test 
site  this  ratio  is  2.  For  an  actual  underground  storage  site  a greater 
ratio  might  be  favoured. 

Table  4.2  shows  the  total  effect  of  these  phenomena  as  estimated  by  the 
model . 

For  a loading  density  of  100  kg  TNT/m^  we  can  see  that  the  pressure 
undergoes  a total  decrease  of  about  50  % while  it  propagates  through 
the  branch  tunnel  and  the  junction  into  the  transport  tunnef.  This 
demonstrates  that  this  connection  between  the  chamber  and  the  tran- 
sport tunnel  constitutes  an  effective  blast  trap. 


5 CONCLUSION 


5.  1 Expected  effects  in  a tunnel  outside  neighbouring  chambers  in  the 

event  of  a mass  detonation  in  a storage  chamber 

In  the  full  scale  case  and  with  a loading  density  of  200  kg  TN'f/m^ 
it  is  concluded  that  the  air  shock  will  probably  be  within  the  range 
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55  - 80  bars.  If  the  loading  density  is  reduced  to  100  kg  TNT/rn^ 
the  air  shock  overpressure  at  the  entrance  to  the  next  chamber  is 
expected  to  be  within  a range  of  35  - 55  bars. 

The  air  shock  wave  will  pass  the  entrance  to  a neighbouring  chamber 
about  15-25  ms  after  the  detonation.  The  explosive  gases  follow 
behind  the  air  shock  at  temperature  expected  to  be  above  1000  °K. 

Blast  waves  following  shots  fired  at  the  end  of  the  tunnel  may  within 
about  10  tunnel  diameters  from  the  charge  have  a peak  pressure  of 
3 to  5 times  the  peak  pressure  in  blast  waves  from  corresponding 
charge  sizes  fired  in  the  chamber. 

The  peak  pressure  in  the  blast  wave  following  a detonation  of  a 1000  kg  TNT 
charge  located  close  to  the  end  of  the  tunnel  is  demonstrated  to  be  con- 
siderably higher  than  the  peak  pressure  in  the  blast  wave  following 
a detonation  of  a similar  charge  placed  a few  tunnel  diameters  from 
the  tunnel's  end.  The  attenuation  of  the  first  mentioned  blast  wave  is, 
however,  greater  in  this  case  so  the  peak  pressure  of  this  blast  wave 
may  with  increasing  distance  decay  to  a value  below  the  peak  pressure 
of  the  latter  at  corresponding  places  in  the  tunnel. 

The  term  duration  is  not  thought  to  be  particularly  useful  to  blast  waves 
in  tunnel  systems.  For  a full  scale  detonation,  however,  the  time  in 
which  an  appreciable  pressure  is  measured,  is  in  the  order  of  1 second. 

Installation  of  blast  doors,  as  described  in  Ref  (6),  at  the  entrance  to 
the  storage  chambers  in  underground  ammunitions  storage  3itcs  similar 
to  the  test  site,  will  probably  provide  effective  protection  against  the 
air  shock  and  the  explosive  ga9es  to  the  chambers  behind  for  loading 
densities  higher  than  100  kg  TNT/m^. 
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Shot  2 

Shot  3 

Shot  4 and  5 

Shot  6 

Dist 

m 

Time 

ms 

Dist 

m 

Time 

m3 

Dist 

m 

Time 

ms 

Dist 

m 

Time 

ms 

IB.  90 

16.  5 

21.15 

21.1 

21.20 

20.0 

13.90 

15.6 

21.10 

29.  5 

31.35 

38.0 

31.40 

35.3 

29.10 

24.0 

39.10 

43.0 

59.90 

93.4 

41  40 

52.1 

39.  10 

36.8 

49. 10 

57.6 

80.42 

140.  0 

51.40 

69. 1 

49.10 

47.0 

59.95 

84.5 

57.65 

87.5 

78.17 

90.  5 

Shot  7 

Shot  8 

Shot  9 

Shot  10 

Dist 

m 

Time 

ms 

Dist 

m 

Time 

ms 

Dist 

m 

Time 

ms 

Dist 

m 

Time 

ms 

21.20 

20.0 

21.85 

15.  5 

17.60 

11.6 

21.90 

18.  9 

31.40 

32.0 

32.05 

22.3 

37.  80 

27.8 

42.10 

36.  5 

41.40 

48.0 

42.05 

30.2 

47,80 

36.0 

52.10 

44.5 

51.40 

62.7 

52.05 

38.  5 

58.3  5 

45.  5 

60.65 

54.0 

59.95 

75.8 

60.60 

47.4 

76.  87 

71.0 

81.17 

80.6  ' 

81.12 

75.0 

Table  3 . 2 Measured  arri  /a!  time  distance  data 


Shot 

2 

3 4 

and  5 

6 

7 

8 

9 

10 

Distance 

Time 

Time 

Time 

Time 

Time 

Time 

Time 

Time 

m 

rns 

ms 

ms 

ms 

ms 

ms 

ms 

ms 

20.00 

17.86 

19.46 

18.14 

16.44 

ie.i3 

1 1.65 

13.27 

17.71 

30.00 

30.63 

35.45 

33.34 

25.59 

31.12 

20.62 

20.63 

25.07 

40.00 

44.33 

53.42 

49.52 

36.26 

45. 18 

28.01 

29.15 

33.52 

50.00 

58.92 

73.03 

66.63 

18.42 

60.25 

36.93 

38.86 

43.11 

60.00 

71.35 

93.96 

34 . 58 

62.01 

76.26 

47.40 

49.76 

53.83 

70.00  , 

90.57 

115.95 

103.30 

76.93 

93.16 

59.41 

61.64 

65.68 

80.00 

1 07 . 53 

133.83 

122.71 

93.09 

110.86 

72.90 

75.06 

78.63 

St.  deviation 

0.38 

0.13 

0.  1 

0.  fci 

0.  50 

0.  SO 

0.46 

0.74 

Function  used 

for  curve 

fitting  : S/C-T  = 

* A + A . 

o 1 

Ln  (T+Tq) 

S = distance, 

T = time. 

C = sound  velocity  at  7 °C 

i Aq  , and  T 

are  constants  to  be  determined 

for  each  fitting 


Table  3.3  Arrival  time  distance  data  obtained  by  curve  fitting  to  the  measured  data  of  Table  3,2 
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Charge 

Charge 

| Initial 

P 

eak  pressure  at 

Peak  pressure  at 

Measured  peak  pressure 

Shot 

weight 

density 

tpre  ssure 

inlet  of  . outlet  of 

t 

.transport 

inlet  of 

outlet  of 

transport 

according 

according 

calcu- 

chamber 

branch 

branch 

tunnel 

branch 

branch 

tunnel 

to  inter- 

to inter- 

lated 

tunnel 

tunnel 

tunnel 

tunnel 

pret  I 

pret  II 

pressure 

kg 

kg/m 

bars 

bars 

bars 

bars 

bars 

bars 

bars 

bars 

bars 

bars 

3 

100 

2 

29 

1C 

9 

6 

10 

9 

6 

5.5 

3.2 

3.4 

4 

300 

7 

62 

16 

11 

10 

15 

13 

9 

7.6 

6.4 

3.6 

5 

300 

7 

62 

16 

14 

10 

15 

13 

9 

8.1 

5.3 

3.6 

7 

1000 

20 

190 

29 

21(25) 

17(16)  | 

26 

21(22) 

15(16) 

1? 

12.5 

5.4 

10 

s;oo 

100 

1070 

69 

50(55) 

36(10) 

58 

43(17) 

31(34) 

32 

20.1 

200 

i 2900 

120 

79(89) 

56(61) 

99 

67(76) 

19(54) 

I 

1 

! 

i 

270 

| 

•1200 

l 

144 

90(103) 
1 

68(72)  J 

117 

1 

' 

76(87) 

i 

55(63) 

1 

! 

rj 

CD 

1 

The  cross-section  area  of  the  branch  tunnel  is  2.4  m . The  numbers  in  the  parenthesis  indicate 

2 

what  to  be  expected  if  the  cross-section  area  of  the  branch  tunnel  is  5 m 


Table  4 . 2 Comparison  of  model  with  measurements 
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Figure  3.3  Measurement  station  situated  in  the  roof 
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Fig  3.5  Blast  pressure  derivel  from  the  arrival  times  for  charges 
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Shot  3 - 100kg  TNT- Chamber 


Fig.  3.10  Measured  overpressure 
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Shot  4 - 300kg  TNT -Chamber 
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Shot  5-  300kg  TNT-Chamber 


Fig.  3.12  Measured  overpressure 
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Shot  5 -300kg  TNT  - Chamber 


Fig.  3 *3  Measured  overpressure 
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Shot  6 - 300kg  TNT -Tunnel 
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Shot  6 - 300kg  TNT -Tunnel 


Fig.  3.15  Measured  overpressure 
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Shot  7 - 1000kg  TNT -Chamber 
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Shot  7 - 1000kg  TNT  - Chamber 
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Fig.  3. 17  Measured  overpressure 
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Fig.  3.21  Measured  overpressure 
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